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Введение. Рассмотрен процесс формирования усталостных 
дефектов в металлических сплавах с различной структур- 
ной морфологией. Цель работы — создание расчетного 
аппарата для определения момента зарождения указанных 
дефектов в условиях циклического нагружения. 
Материалы и методы. Построена физическая модель, 
представлены расчетные выражения. В основу физической 
модели положена теория дислокаций. Показано, что при 
динамических циклических нагружениях определяющее 
значение в процессе зарождения очагов разрушения имеет 
структурный фактор. В зависимости от структуры и 
свойств материала, а также от характера нагрузок, крити- 
ческий усталостный дефект развивается в форме трещины, 
поры или микрократера износа. 

Результаты исследования. Выполнен численный экспе- 
римент по определению момента зарождения дефекта 
критического размера в сплавах на основе железа при 
высокоскоростных капельных соударениях. Представлены 
сравнительные данные расчетов и стендовых испытаний 
по каплеударной эрозии сталей и сплавов со структурой 
феррита, аустенита, сорбита и мартенсита. Оценен вклад 
стадии зарождения в инкубационный период эрозионного 
износа исследованных материалов. 

Обсуждение и заключение. Строгие инструментальные 
методы для определения продолжительности стадии за- 
рождения отсутствуют, поэтому рекомендовано использо- 
вать предложенную расчетно-аналитическую модель. 
Кроме того, выполненная работа дала важный прикладной 
результат — продемонстрировала, что целенаправленное 
конструирование структуры материала может существен- 
но увеличить износостойкость. 


Ключевые слова: сплавы на основе железа, субструктура 
сплавов, циклические нагрузки, зарождение очагов раз- 
рушения, физическая модель, численный эксперимент, 
каплеударная эрозия. 
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Введение. Моделирование специфических непрогнозируемых явлений весьма актуально с точки зре- 
ния предотвращения их нежелательных последствий. Например, в аэронавтике, теплоэнергетике, газотурбин- 
ном производстве и сельском хозяйстве остается нерешаемой проблема каплеударной эрозии. Вопрос обсужда- 
ется на мировом уровне [1-4]. Первые успехи в моделировании этого явления касались единичных капельных 
соударений [5—10]. В некоторых работах [11-13] предприняты попытки количественно оценить способность 
материалов и покрытий сопротивляться действию динамических циклических нагрузок. В частности, удалось 
создать расчетно-аналитическую модель, апробированную в условиях жидкокапельного ударного воздействия 
[14-18]. Многие авторы [19-21] не дифференцируют поведение материала или покрытия на стадиях зарожде- 
ния и развития разрушения. Зарождение весьма проблематично идентифицируется в эксперименте. В названной 
выше модели оно оценивалось с помощью эмпирических коэффициентов, не имевших универсального теорети- 
ческого обосновании, и это явный недостаток решения. Цель настоящей работы — более глубокое теоретиче- 
ское рассмотрение стадии зарождения дефектов в покрытии, приводящих к его износу и разрушению, а также 
создание на этой основе расчетного аппарата для определения момента зарождения трещин или пор в различ- 
ных материалах и покрытиях в условиях циклического нагружения. 

Материалы и методы. Трудности оценки стадии зарождения разрушения и износа при циклических 
воздействиях (например, при высокоскоростных жидкокапельных соударениях, провоцирующих каплеударную 
эрозию металлических изделий) связаны не только с инструментальным определением зародышей трещин или 
пор. На этой стадии наличие структурных элементов в материале не позволяет представить процесс непрерыв- 
ным континуумом и использовать для расчетов линейную механику разрушения (например, усталостную тео- 
рию Пэриса — Эрдогана). Поэтому в механике разрушения начальная стадия зарождения усталостных дефек- 
тов выделяется в особый этап. Он называется «микроструктурная механика разрушения» и протекает вплоть до 
достижения дефектом размера /,= (4...10)-4, где 4 — размер структурного элемента материала [22]. Аналитиче- 
ская оценка стадии зарождения усталостных дефектов требует, прежде всего, металлофизического подхода. 
Задача заключается в нахождении расчетного выражения для числа циклических нагружений №, необходимых 
для формирования в материале усталостного дефекта критического размера. В качестве эмпирической базы 
расчетной модели можно использовать также экспериментальные данные стендовых испытаний на каплеудар- 
ную эрозию, где роль № выполняет количество капельных соударений. 

В общем случае уравнение для числа динамических циклических нагружений №; должно включать три 
составляющие: механическую, кинетическую и структурную. 

Металлофизическая концепция механической составляющей связана с числом (плотностью) подвиж- 
ных дислокаций р», возникающих при единичном ударном нагружении (соударении). В процессе многократ- 
ных соударений количество подвижных дислокаций увеличивается, они перемещаются в металлической мат- 
рице по плоскостям скольжения под действием касательных напряжений с, до тех пор, пока не сформируют на 
ближайших непреодолимых барьерах плоские скопления с критической плотностью ру,. Превышение ру. ведет 
к самопроизвольному разрыву межатомных связей в металле и образованию зародыша трещины. Значения ри. 
для различных материалов и покрытий известны [23]. Механическая составляющая числа соударений выража- 
ется в форме зависимости № = / (ру, м). 

Кинетическая составляющая числа соударений № учитывает зависимость энергии дислокации от ее 
скорости Г.. Из теории дислокаций известно, что по мере возрастания скорости дислокации ее энергия увели- 
чивается в соответствии с выражением Эйнштейна для тел, движущихся со скоростями, близкими к скорости 
света. Только пределом для скорости дислокации является скорость звука в кристалле Со, при которой энергия 
дислокации становится бесконечной. Поэтому с учетом специфики каплеударного воздействия кинетическая 
составляющая числа соударений выразится как № = 1 (Тл/ Со). 

Структурная составляющая числа соударений № имеет два аспекта. Первый заключается в том, что 
движение дислокаций при ударе ограничено размерами структурного элемента матрицы ДО, в пределах которого 
возможен свободный пробег дислокаций. В общем случае в качестве такого элемента рассматривается размер 
зерна, а для тонких покрытий — его толщина. Второй аспект учитывает, что в пределах зерна или покрытия 
могут существовать препятствия для движения дислокаций: частицы вторичных фаз, малоугловые границы, 
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закрепленные атмосферами неподвижные дислокации, дефекты упаковки, сопротивление решетки (напряжение 
Пайерлса). Как отмечено выше, оба эти аспекта отражены в структурной составляющей числа соударений, и 
она представляет собой функцию двух переменных № = (О, АС), где АС — свободная энергия Гиббса для ак- 
тивации процесса преодоления препятствий, находящихся в пределах структурного элемента О. Механизм пре- 
одоления зависит от природы препятствия, что и определяет значение ДС для каждого конкретного случая 
структуры материала (их классификация приведена в [23]). 

Из сказанного следует, что все составляющие искомого числа соударений № влияют на процесс фор- 
мирования усталостного дефекта одновременно — при каждом соударении. Этот факт определяет коммутатив- 
ный характер взаимодействия механической, кинетической и структурной составляющих. Тогда в наиболее 
общем виде № определяется выражением 


ие р — 
ма 1 [2]. .е АТ, (1) 





Рт Со № 
где А — постоянная Больцмана; Т — термодинамическая температура, К; /› — путь, пройденный подвижными 
дислокациями за один цикл нагружения (соударение). 

Выражение (1) получено на основе теории дислокаций и представляет собой теоретическую концеп- 
цию. Для использования (1) при оценке износостойкости материалов и покрытий в условиях циклических жид- 
кокапельных соударений необходимо раскрыть входящие в него величины (такие как ри, Га, 1, АС) через изме- 
ряемые параметры соударения (например, скорость соударения Го и размер капли Ао). Авторы уже выполнили 
соответствующую детализацию модели (1) и готовят материал к публикации. 

Результаты исследования. Расчетная модель апробирована для сплавов на основе железа различного 
состава и структуры. Такой выбор обусловлен, прежде всего, различием в поведении дислокаций в этих сплавах 
и позволяет выявить возможности расчетной модели. Так, в феррите и аустените фактор, тормозящий дислока- 
ции, — это сопротивление решетки (напряжение Пайерлса). В сорбитной структуре скольжение дислокаций 
определяется морфологией и распределением препятствий — карбидами Ее и Сг. В мартенсите движение дис- 
локаций затруднено изначально. Таким образом, предметом исследования этой части работы фактически явля- 
ется структурный фактор. 

Поскольку все исследованные материалы представляют собой сплавы на основе железа, их основные 
физико-механические характеристики отличаются не существенно. В расчетах использовались следующие зна- 
чения: модуль упругости Ё= 186...218 ГПа; коэффициент Пуассона у= 0,20...0,31; модуль сдвига 
ц= 64...80 ГПа; вектор Бюргерса Ь = 2,5.10`'° м. Данные численных экспериментов по реализации представ- 
ленной модели и стендовых испытаний сплавов на основе Ее с различной структурной морфологией, в том чис- 
ле и гетерогенных, приведены в табл. | и на рис. 1. 





























Таблица 1 
Экспериментальные т и расчетные №; значения числа соударений* для сплавов на основе Ее 
Параметры соударения 
Параметры структуры 
Материал Г, = 250 м/с; В, = 0,55 мм Г, = 340 м/с: Юу= 0,32 мм 
структура Р Расчет Эксп. Расчет Эксп. 
(структура) ешетка аа |ры$| м8] 
(М,) № 0% то № 0 то 
АРМКОЕ 
р ОЦК (2,9) 0,5 100 3,0 2197 — — 1187 — — 
(феррит) 
08Х18Н10Т 
ГЦК (3,06) | 0,5 100 3,0 2446 0,253 9680 1322 0,154 8597 
(аустенит) 
20х13 0,35— 
К (2 2 1 | 4 2 1 42 254 14 
еж ОЦК (2,9) ‚0 00 0,50 655 0,280 6630 37 0,25 768 
20х1 Т | 
О а | 0,02 | 0,01 | 0,01 | 9844 | 0,550 | 17861 | 5705 | 0,470 | 12140 
(мартенсит) (2,95) 



































*Представлены значения чисел капельных соударений, попадающих в одну точку (определяются в программ- 





ном обеспечении испытательного стенда). 








Кудряков О. В. и др. Анализ начальной стадии усталостного износа гетероструктурных материалов в условиях 


КиагуаКоу О. У., ег а1. Апа[у515 ор йе тии 5азе ор йвие теаг т пеегоягисиге тщена!5 ип4ег сотасй сусПс оа@ тя 

















ыы Ей 
с 












































ы 
8 
— 





10 20 30 40 
И,, тб 7, т 


а) 6) 






































10 20 30 40 50 60 10 
Г, тб 


в) 
Рис. 1. Зависимости числа соударений №, необходимых для зарождения трещин, от скорости соударения 
И, при Ко=0,55 мм для сплавов на основе Ее с различной структурой: феррит железа АРМКО (а); сорбит стали 20Х13 (6); 
мартенсит стали 20Х13 (в) 


Значения некоторых параметров, приведенных в табл. 1, нуждаются в пояснении. В качестве расстоя- 
ния между препятствиями [ на пути подвижных дислокаций В феррите и аустените принимается расстояние 


между дислокациями, то есть [= у Ра . Учитывая высокую пластичность этих твердых растворов, плотность 


их дислокаций относительно невелика р. = 10" м’, что дает [= 3 мкм. В закаленной стали плотность дислока- 
ций близка к критической р», = 10° м”, что не только дает значение структурного параметра /= 0,01 мкм в 
табл. 1, но и изменяет методику расчета числа соударений в исходной модели (1). 

Для упрочненных металлических материалов (например, для сильно деформированных или закаленных 
сплавов, в том числе 20Х13 с мартенситной структурой), в выражении (1) первые три сомножителя обращаются 
в единицу. В таких материалах критическая плотность дислокаций р/у. уже достигнута, поэтому понятие по- 


движных дислокаций ри теряет смысл. В численном виде это выражается так: р„„/р„=1; 7.=0; Р=Ц = у х [р Е 


Что же касается экспоненты в выражении (1), то она сохраняет свой смысл со значениями параметров 0 и ГА 
приведенными в табл. 1. Чтобы в такой структуре зародилась трещина, внешнему циклическому воздействию 
необходимо преодолеть существующие напряжения в насыщенной дислокационной среде и сформировать де- 
фект сдвигового характера. Учитывая сказанное, расчетное выражение (1) ДЛЯ закаленной стали примет 


вид [23—26]: 
№; = Се ы т. (2) 


0; 

Здесь АГ — энергия активации процесса преодоления препятствия без приложения внешнего напряжения, ве- 
личина ДЕ фактически определяет прочность препятствий с точки зрения прохождения их дислокацией; о. — 
внутреннее напряжение, существующее в материале и позволяющее дислокации проходить препятствие при 
минимальных значениях ДС. В качестве о. принимается напряжение течения твердого тела при Т = ОК. Вели- 
чины о. и ДЕ являются свойствами материала и в общем случае выражаются через его базовые физические ха- 
рактеристики ци В: 

о. = И.В /ЕАР=е. и. В (3) 
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Коэффициент о в выражениях (3) классифицирует препятствия по их прочности [23]. 

В качестве напряжения с; в выражении (2) рассматривается сдвиговое напряжение в плоскости сколь- 
жения, под действием которого происходит перемещение подвижных дислокаций. Напряжение с, представляет 
собой проекцию вектора внешнего напряжения на плоскость скольжения и связано с системой скольжения ре- 
шетки через фактор Тейлора М, по закону Шмидта — Боаса для поликристалла [27]. Оно может быть выражено 
через параметры соударения Гу и Ки [28]. 

Экспериментальные значения, показанные в табл. 1, характеризуют инкубационный период эрозионно- 
го износа материала то. Эта величина включает как стадию зарождения дефекта (трещины, поры, микрократе- 
ры), так и стадию его развития до начала потери массы образцом. Не представляется возможным инструмен- 
тально выделить в экспериментальных величинах ту количество соударений, приходящихся только на стадию 
зарождения, поэтому значения № являются исключительно расчетными. Вклад стадии зарождения в общее зна- 
чение характеризуется коэффициентом 04 = № / то. 

Кроме численных значений, представленных в табл. 1, модель учитывает количество соударений №, 
необходимое для зарождения очага разрушения (трещины, поры, микрократеры), и скорость соударений Г, при 
фиксированных значениях размера капли Ко. Кроме того, используя данное решение, можно определить зави- 
симость № от Ио. Некоторые варианты таких зависимостей представлены на рис. 1. 

Полученные результаты показывают, что расчетные данные Л; не противоречат эксперименту (данным 
стендовых испытаний 770). Они также демонстрируют соответствие базовым канонам механики разрушения. В 
частности, коэффициент 0% характеризует энергию зарождения по отношению к энергии разрушения. В пла- 
стичных металлических материалах (АРМКО, аустенит, сорбит) энергия зарождения усталостных трещин су- 
щественно меньше энергии их роста. В упрочненных сплавах (мартенсит) энергия зарождения практически все- 
гда выше энергии роста. Как правило, это соотношение значительно больше в пользу стадии зарождения. Дан- 
ные табл. 1 показывают, что в стали 20Х13 с мартенситной структурой стадия зарождения усталостного дефек- 
та № составляет около половины всего инкубационного периода формирования износа ти: 0% = 0,47 и 0,55 для 
скоростей соударения Г, = 340 м/с и 250 м/с соответственно. Тогда как для той же стали с сорбитной структу- 
рой коэффициент 0% вдвое меньше. Полученные нами значения од соответствуют данным других авторов для 
аналогичных условий циклического нагружения (например, для стали ЗОХГСН2А в [29]. 

Расчетные зависимости на рис. | также соответствуют данным экспериментов. Они расположены 
асимптотически по отношению к значениям И = 100...150 м/с (проведенный численный эксперимент эту об- 
ласть подробно не исследовал). Названные выше расчетные зависимости соответствуют итогам стендовых ис- 
пытаний, которые показали, что при Ио < 135 м/с в сплавах на основе железа износ вообще не наблюдается [30-— 
34]. То есть асимптотичность графиков на рис. 1 подтверждает этот эмпирический факт: при указанных низких 
скоростях соударений не образуются усталостные дефекты критического размера. 

Обсуждение и заключение. Не существует инструментальных методов точного определения продол- 
жительности стадии зарождения очагов разрушения материалов, поэтому при необходимости используется 
предложенная расчетно-аналитическая модель. Она применима к материалам различной морфологии, в кото- 
рых механизмы зарождения дефектов (трещин, пор, микрократеров) имеют дислокационную природу. Поэтому 
при нормальных условиях модель не применима к керамическим материалам с высокой долей ковалентных 
связей. Другим ограничением использования модели является, пожалуй, слишком высокая чувствительность 
экспоненциального множителя в выражениях (1) и (2). Небольшая погрешность при определении численных 
значений величин 041, о, или / в этих формулах не позволяет получить результаты, адекватные эксперименталь- 
ным данным. Выполненная работа дала важный прикладной результат — продемонстрировала, что целена- 
правленное конструирование структуры материала может существенно увеличить износостойкость. 
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